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ABSTRACT 

We present an analysis of the diffuse X-ray emission in 19 compact groups (CGs) of galaxies observed with 
Chandra. The hottest, most X-ray luminous CGs agree well with the galaxy cluster X-ray scaling relations in L x -T 
and L x -a , even in CGs where the hot gas is associated with only the brightest galaxy. Using Spitzer photometry, 
we compute stellar masses and classify Hickson CGs 19, 22, 40, and 42, and RSCGs 32, 44, and 86 as fossil 
groups using a new definition for fossil systems that includes a broader range of masses. We find that CGs with 
total stellar and Hi masses >10 11,3 M© are often X-ray luminous, while lower-mass CGs only sometimes exhibit 
faint, localized X-ray emission. Additionally, we compare the diffuse X-ray luminosity against both the total UV 
and 24 /zm star formation rates of each CG and optical colors of the most massive galaxy in each of the CGs. The 
most X-ray luminous CGs have the lowest star formation rates, likely because there is no cold gas available for 
star formation, either because the majority of the baryons in these CGs are in stars or the X-ray halo, or due to 
gas stripping from the galaxies in CGs with hot halos. Finally, the optical colors that trace recent star formation 
histories of the most massive group galaxies do not correlate with the X-ray luminosities of the CGs, indicating 
that perhaps the current state of the X-ray halos is independent of the recent history of stellar mass assembly in the 
most massive galaxies. 
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1. INTRODUCTION 

In the local universe, the majority of galaxies exist in gravi- 
tationally bound systems, i.e., in groups or clusters (e.g., Tully 

1987; Small et al. 1999; Karachentsev 2005). Cosmological 
ACDM models would imply that in rich groups and clusters of 
galaxies the fraction of baryons in stars may be as little as 20% 
(Borgani et al. 2004), while observations indicate a value closer 
to 10% (e.g., Balogh et al. 2001 ; Lin et al. 2003). The remaining 
baryons are in the form of gas in various states, i.e., molecular, 
neutral, or ionized. In rich clusters, ram-pressure stripping and 
harassment of gas-rich galaxies deposits vast quantities of neu- 
tral gas into the intracluster medium (ICM), providing material 
for a virialized X-ray halo (cf., Gunn & Gott 1972). 

In galaxy clusters, the hot ICM is already largely developed, 
therefore we must look to the building blocks of clusters to 
examine the early stages of the growth of the hot gas halos. 
Compact groups (CGs) have high galaxy number densities 
similar to the cores of rich clusters, and they are expected to 
experience enhanced tidal encounters and mergers compared 
to loose groups while their low velocity dispersions lengthen 
the timescales over which these encounters occur relative to 
clusters. These systems provide excellent laboratories to study 
the effects of galaxy interactions on the build-up of hot gas halos 
in low-mass groups of galaxies, which are the building blocks 
of rich clusters (Peebles 1970; Gonzalez et al. 2005). 

Several studies have investigated the X-ray properties of CGs 
and specifically the diffuse, hot gas in these systems (e.g., 
Helsdon et al. 2001; Desjardins et al. 2013; Fuse & Broming 


2013). The first comprehensive examination of diffuse X-ray 
emission in CGs was performed by Ponman et al. (1996), 
in which the authors use X-ray observations with the ROSAT 
Position Sensitive Proportional Counter (PSPC) to examine the 
group-linked hot gas in a sample of 85 Hickson CGs (HCGs; 
Hickson 1982) of which 22 were detected. The authors made 
efforts to mask the soft X-ray emission from the individual 
galaxies and report that the remaining emission appears to be 
clumpy, suggesting that, in contrast to clusters, the hot gas is 
not in equilibrium. 

With the much improved spatial and spectral resolution of the 
Chandra Advanced CCD Imaging Spectrograph (ACIS) com- 
pared to the ROSAT PSPC (angular resolutions 075 and 25" 
FWHM, respectively), Desjardins et al. (2013) find that the de- 
tectable diffuse X-ray emission in a small sample of nine HCGs 
have varied morphologies that range broadly from linked to the 
individual galaxies to a true intragroup medium (IGM; not to 
be confused with the intergalactic medium). The galaxy-linked 
emission is typically associated with vigorous star formation, 
while hot gas in the form of an IGM is likely due to virialization 
of the baryons by the group potential well. HCG 42 may be an 
exception as it has a hot gas halo associated with the brightest 
group galaxy resembling a hot IGM that appears small in ex- 
tent, but which may extend farther than is detectable due to low 
surface brightness. 

In this study, we expand upon the analysis presented by 
Desjardins et al. (2013) using Chandra ACIS archival observa- 
tions of an additional 10 CGs, thus bringing the total Chandra 
sample to 19 CGs when combined with Desjardins et al. (2013). 
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Table 1 

Compact Group Sample 


Desjardins et al. 


Group 

Coordinates (J2000) 

z 

VCMB a 

<T b 

R c 

M H, 

Afdyn 

References 


a 

8 


(km s -1 ) 

(km s -1 ) 

(kpc) 

(10 9 Mo) 

(10 11 Mo) 


HCG 30 

04 h 33 m 28 s 

— 02°49'57" 

0.0154 

4562 

121-23 

81.2 ±2.2 

0.60 ± 0.06 

13.80 ±0.37 

1,3, 4, 7 

HCG 37 

09 h 13 m 35 s 

+30°00'51" 

0.0223 

6940 

451 ± 17 

42.4 ± 2.8 

5.40 ±0.15 

100.257^ 

1,3,4, 8,9 

HCG 40 

09 h 38 m 54 s 

— 04°51'07" 

0.0223 

7026 

148 ± 15 

37.6 ± 0.8 

6.60 ±0.15 

9.63 ± 0.21 

1,3,4, 10, 11 

HCG 51 

1 l h 22 m 21 s 

+24°17'35" 

0.0258 

8051 

26871 

93.7 ± 3.2 

<4.57 

78.39 ± 2.68 

1,3,5, 8, 12 

HCG 68 

13 h 53 m 41 s 

+40°19'07" 

0.0080 

2583 

10871 

60.0 ± 2.2 

5.60 ± 0.02 

7.21 ± 0.30 

1,2, 3,4, 13, 14, 15, 16 

HCG 79 

15 h 59 m 12 s 

+20°45'31" 

0.0145 

4439 

2657 12 

16.5 ± 1.0 

4.20 ±0.13 

13.41 ±0.81 

1,2, 3,4, 10, 17, 18 

HCG 97 

23 h 47 m 23 s 

— 02°19'34" 

0.0218 

6174 

359 ± 18 

99.8 ± 5.9 

4.3 ±0.13 

149.637I 87 

1,2, 3,4,9, 11, 19 

HCG 100 

00 h 01 m 21 s 

+13°07'57" 

0.0178 

4976 

142 ± 24 

54.1 ± 1.5 

9.0 ± 0.09 

12.65 ± 0.36 

1,3,4, 20,21 

RSCG 17 

01 h 56 m 22 s 

+05°38'37" 

0.0190 

5415 

386 ± 19 

35.3 ± 2.4 

<15.85 

61.147 4 4 2 2 3 1 

2, 6, 9, 12 

RSCG 31 

09 h 17 m 22 s 

+41°57'24" 

0.0060 

2009 

52 ± 12 

41.3 ± 1.1 

1.82 

1.28 ±0.04 

2, 6, 22, 23, 24 


Notes. 

a Value taken from the NASA Extragalactic Database (NED) based on Fixsen et al. (1996). 
b We recalculated the velocity dispersion o using the CG member velocities from the cited references. 
c The median projected two-galaxy separation. 

References. Catalogs: (1) Hickson 1982; (2) Barton et al. 1996 — positions; (3) Hickson et al. 1992 — Hi masses: (4) Borthakur et al. 2010; (5) Verdes-Montenegro 
et al. 2001; (6) L. M. Walker et al. (in preparation) — velocities: (7) Jones et al. 2009; (8) Falco et al. 1999; (9) Mahdavi & Geller 2004; (10) Nishiura et al. 2000; (11) 
de Carvalho et al. 1997; (12) de Vaucouleurs et al. 1991; (13) van Driel et al. 2001; (14) Denicolo et al. 2005; (15) Rhee & van Albada 1996; (16) Simien & Prugniel 
2002; (17) Paturel et al. 2003; (18) Bonfanti et al. 1999; (19) Mahdavi et al. 2005; (20) Huchra et al. 1999; (21) Hickson 1993; (22) Strauss et al. 1992; (23) Monnier 
Ragaigne et al. 2003; (24) Nordgren et al. 1997. 


We also utilize Spitzer , Sloan Digital Sky Survey (SDSS), and 
Apache Point Observatory (APO) data to characterize how the 
total CG stellar mass relates to the observed X-ray properties 
of the groups. Section 2 describes the sample selection and the 
general characteristics of the CGs in our study. In Section 3, we 
list the Chandra data reduction steps and determine group stellar 
masses. Section 4 discusses our findings and their implications, 
and Section 5 summarizes our conclusions. Errors are reported at 
the 90% confidence level unless otherwise stated. For all calcu- 
lations, we assumed the currently favored cosmological param- 
eters of Q m = 0.27, Qa = 0.73, and Ho = 70 km s -1 Mpc -1 
(Hinshaw et al. 2013). 

2. SAMPLE DESCRIPTION 

The selection criteria of HCGs and Redshift Survey CGs 
(RSCGs; Barton et al. 1996) yield samples that are well 
suited for our study of the hot gas properties in dense galaxy 
groups. The HCG catalog was compiled by searching Palomar 
Observatory Sky Survey (POSS) data, which cover the sky at 
all declinations 8 > —27°, for groups with N > 4 galaxies 7 
with magnitudes within 3 mag of the brightest group galaxy, 
On ^ 3 6g, and \±g < 26 mag arcsec -2 . In this definition, 
0 N is the angular diameter of the largest concentric circle 
that does not include non-member galaxies within 3 mag of 
the brightest group galaxy, 0q is the angular diameter of the 
smallest concentric circle that includes the nuclei of all of the 
group members, and \±g is the surface brightness averaged over 
the circle defined by 0q . The last two criteria for the HCG 
catalog pertain to the isolation and compactness of the groups, 
respectively. All photometric measurements were made in the 
POSS E band (most equivalent to the standard Johnson R filter), 
for which the POSS observations are complete to m = 20.0 mag. 

Barton et al. (1996) used the second Center for Astrophysics 
redshift survey (CfA2) and second Southern Sky Redshift 


7 This requirement has been relaxed to N ^ 3 due to the discovery as a result 
of spectroscopic follow-up that only 69% of HCGs have N ^ 4 accordant 
members, while 92% have N ^ 3 (Hickson et al. 1992). 


Survey (SSRS2) to identify a sample of CGs with similar 
properties to the Hickson (1982) catalog. The CfA2 consists of 
a strip that covers approximately 1 17° x 6° and is centered near 
the north Galactic pole (e.g., de Lapparent et al. 1986), while 
the SSRS2 covers 1.3 sr around the southern Galactic cap (da 
Costa et al. 1994). Both surveys are complete tomB 0 = 15.5 and 
Barton et al. (1996) only consider galaxies in the line-of-sight 
velocity range 300 < v < 15, 000 km s _1 . These RSCGs are 
selected using a friends-of-friends algorithm in which groups 
with N ^ 3 are considered RSCGs, and group members are 
found using Vo ^ 1000 km s -1 and Do ^ 50 kpc, where Vo and 
Do represent the velocity difference and projected separation, 
respectively, between neighbor galaxies. Barton et al. (1996) 
chose the value of Do to most closely match the observed 
properties of the HCG sample. 

Our sample consists of all HCGs and RSCGs available in the 
Chandra data archive that are not part of our previous diffuse 
X-ray study presented in Desjardins et al. (2013) and that are 
completely covered by the Chandra footprint. We then refine our 
sample by comparing RSCGs with nearby galaxy clusters and 
removing those groups with projected separations of <1 Mpc 
and velocity differences <3cr from a galaxy cluster. The group 
mean position and velocity were used for comparison. This 
is necessary because we wish to include only isolated CGs, 
and the RSCG catalog does not employ an isolation criterion 
similar to that of the HCGs. Due to the lack of such a criterion, 
the RSCG catalog includes a number of dense environments 
misidentified as CGs, e.g., RSCGs 67 and 68 are in the core 
of the Coma Cluster. The criteria described above resulted 
in a sample of 10 additional CGs in our study compared to 
Desjardins et al. (2013), which we list in Table 1. We note, 
however, that Desjardins et al. (2013) did include HCG 42, 
which is a high galaxy density region potentially located in a 
filamentary structure along the line of sight (Konstantopoulos 
et al. 2013). Indeed, Diaz-Gimenez & Mamon (2010) suggest 
that at best ~85% of HCGs with more than four members are 
truly dense systems rather than chance alignments. This leads 
to an estimate that two to three of the groups in our sample may 
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Table 2 

Extended Group Membership 


Group 

C500 a 

(arcmin) 

Galaxy 

Magnitude 

(mag) 

Ang. Sep. 
(arcmin) 

Reference 

Dwarfs 

HCG 30 

11.0 





0 

HCG 31 

12.5 





2 b 

HCG 37 

11.7 





0 

HCG 40 

7.4 





2 C 

HCG 51 

9.3 





7 

HCG 68 

13.6 





2 

HCG 79 

11.3 





1 

HCG 92 

14.5 





3 

HCG 100 

10.1 

MRK 935 

14.28 (R) 

1.8 

1 

0 

RSCG 31 

25.0 





1 


Notes. 

a Or the angular separation corresponding to a projected 200 kpc radius, 
whichever is larger. 

b One of the two dwarf galaxies in HCG 31 is actually low-mass tidal debris 
designated HCG 31R (Verdes-Montenegro et al. 2005). 
c One of the two dwarf galaxies in HCG 40 (HCG 40-06 in de Carvalho et al. 
1997) has only a 5-band magnitude listed in NED, corresponding to an absolute 
magnitude Mb ^ — 17 mag (de Carvalho et al. 1997), but no 5-band data is 
listed. Based on an inspection of the image, and the fact that the J-, H-, and 
5-band magnitudes from the 2MASS catalog (Skrutskie et al. 2006) are all 
>3 mag fainter than HCG 40A, we classify this as a dwarf group member. 
References. (1) Hickson et al. 1989. 


not be spatially dense, however, some of our systems only have 
three galaxies, which complicates this estimate. Further, our 
sample is not completely random as 16 of the groups between 
our study and Desjardins et al. (2013) were selected to study 
X-ray emission in the group environment (with HCGs 5 1 and 
97 and RSCG 17 chosen because they were known to be X- 
ray bright), while the remaining groups were observed to study 
either X-ray emission in and around early-type galaxies or the 
supernova SN2006jc in RSCG 31. 

Another concern when studying CGs is the possibility of 
additional galaxies far from the compact core, but still bound 
to the group, de Carvalho et al. (1997), Zabludoff & Mulchaey 
(1998), and Konstantopoulos et al. (2010, 2012, 2013) examined 
the extended populations of HCGs 7, 16, 22, 40, 42, 59, 62, 90, 
and 97 and RSCG 17 and found that only HCGs 42, 62, and 90 
and RSCG 17 had substantial populations of galaxies outside 
of the core region, while HCG 97 is missing two relatively 
bright galaxies in the HCG catalog. We therefore label these 
groups as lower-limits with respect to their total group stellar 
masses in the figures throughout this work. For the remaining 
groups in our sample, we used the NASA Extragalactic Database 
(NED) to search for additional galaxies within the larger of r$ oo 
or 200 kpc in radius and ±1000 km s -1 of the group mean 
velocity. The results of this search are show in Table 2. Note that 
we only provide detailed information on the luminous galaxies, 
and simply list the aggregate number of dwarfs in each group. 
To distinguish between luminous and dwarf group members, we 
applied an absolute magnitude cut at — 17 mag in the R band. We 
further required that the galaxy be within three magnitudes of the 
brightest group galaxy to satisfy the Hickson (1982) selection 
criteria, as any galaxies failing this test would likely be of 
little relative importance in determining total group properties. 
Finally, in some cases, photometric data were missing in NED, 
and we examined the images by eye to compare the relative sizes 
of the galaxies on the sky (this resulted in only dwarf galaxy 
classifications). Only HCG 100 excludes a relatively massive 


galaxy in the Hickson catalog, while HCG 5 1 is missing seven 
dwarf galaxies. We label these two additional CGs as lower- 
limits in total group stellar mass in our figures. Note that for 
CGs embedded within larger structures, we only consider the 
properties of galaxies that make up the compact region. While 
this exclusion of the extended populations may seem in error, 
Palumbo et al. (1995) examined the extended populations of the 
Hickson (1982) sample and found that the compact cores and 
extended halos showed statistically different properties (e.g., 
spiral fraction) indicating that the CGs are “disconnected” from 
their environments. Evidence of this distinction between CG 
galaxies and their surrounding environment can be seen in 
the work of Johnson et al. (2007), Walker et al. (2010), and 
Walker et al. (2012) who found a gap in the mid-IR color 
distribution of CG galaxies suggestive of accelerated evolution 
attributed to the CG environment. Further, the galaxies far from 
the compact cores are, in most cases, dwarf galaxies that do 
not add significant stellar mass to the group. Dozens of such 
galaxies would be required to significantly affect our results. 
While the group members far outside the core may also add 
substantially to the total group star formation rate (SFR), these 
members are not yet impacted by ram-pressure stripping nor 
have they contributed much gas to the formation of the IGM, 
therefore we exclude them in the discussion of the link between 
star formation and diffuse X-ray luminosity. 

We also include a comparison sample of galaxy clusters 
from Wu et al. (1999) and Zhang et al. (2011). We se- 
lected clusters from Wu et al. (1999) and Zhang et al. (2011) 
for each of the L x -T , L x -a , and cr-T relationships by 
including only clusters that had published uncertainties for 
both values in each scaling relation (for further details see 
Desjardins etal. 2013). The Wuetal. (1999) clusters are amassed 
from the literature (see their Table 1 for the full list of refer- 
ences), and have redshifts z < 1 and (z) ~ 0.1, temperatures 
1 < T < 17 keV, velocity dispersions 150 < a < 2000 km s -1 , 
and X-ray luminosities 42 < log 10 (Lx) < 46. The Zhang et al. 
(2011) measurements use XMM-Newton observations of 62 
of the 64 Highest X-ray FLUx Galaxy Cluster Sample (“HI- 
FLUGCS”; Reiprich & Bohringer 2002) galaxy clusters, which 
were originally identified using ROSAT X-ray data at Galactic 
latitudes of \t\ > 20°. The clusters from Zhang et al. (2011) 
have luminosities 42 < log 10 (±x) ^$45 erg s -1 , temperatures 
0.7 < T <15 keV, velocity dispersions 200 < o < 1000 km s _1 , 
and a mean redshift of (z) = 0.05. The X-ray properties of 
the Zhang et al. (2011) clusters are measured within r^oo, while 
the Wu et al. (1999) clusters are taken from the literature and 
corrected to a common radius using a ft model. Though it is 
unclear what radius Wu et al. (1999) used, it is reasonable to 
assume this correction was performed to oo- 

3. ANALYSIS 
3.1. X-Ray Observations 

The Chandra observations are summarized in Table 3. Data 
were taken in either FAINT or VFAINT mode with no gratings. 
We performed the calibration of the Chandra data using the 
Chandra Interactive Analysis of Observations (CIAO version 
4.5) tool in conjunction with the CIAO calibration database 
version 4.5.5. 1. Beginning with the Level 1 events file, we 
processed the data using acis_process_events with correc- 
tions for the charge transfer inefficiency and time-dependent 
gain. We used the status bits in the Level 1 events file set 
by the standard data processing pipeline tasks de streak and 
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Table 3 

Summary of Chandra ACIS Observations 


Group 

ObsID Array Mode Exposure 
(ks) 

Date 

Previous 

Publications 

HCG 30 

6977 

S 

VF 

29.7 

2006 Feb 7 

1,2 

HCG 37 

5789 

S 

VF 

17.9 

2005 Jan 13 

1,2 

HCG 40 

5788 

S 

VF 

33.2 

2005 Jan 29 

1,2,3 


6203 

S 

VF 

15.0 

2005 Jan 29 

1,3 

HCG 51 

4989 

S 

VF 

38.5 

2004 Feb 15 

2, 4-10 


5304 

S 

VF 

13.0 

2005 Feb 24 

5, 6,9 

HCG 68 

5903 

S 

VF 

4.5 

2005 Apr 10 

3, 11 

HCG 79 

11261 

S 

VF 

69.2 

2010 May 20 

12 

HCG 97 

4988 

S 

VF 

57.4 

2005 Jan 14 

1,2, 9, 13 

HCG 100 

6978 

I 

VF 

27.8 

2006 Dec 6 

2 


8491 

I 

VF 

17.8 

2007 Jan 24 

1 

RSCG 17 

2223 

s 

F 

30.4 

2001 Jan 28 

4, 5, 7, 9, 10, 14-28 

RSCG31 

6729 

s 

VF 

54.5 

2007 Jan 6 

29, 30 


8457 

s 

F 

9.8 

2006 Nov 4 

11,29-31 


9093 

s 

VF 

24.8 

2008 Jan 20 

30 


10567 

s 

F 

5.1 

2009 Jan 24 

30-32 


References. (1) Rasmussen et al. 2008; (2) Fuse & Broming 2013; (3) Cluver 
et al. 2013; (4) Kim et al. 2007b; (5) Sun et al. 2009; (6) Sun 2009; (7) Dong 
et al. 2010; (8) Haggard et al. 2010; (9) Sun 2012; (10) Trichas et al. 2012; 
(11) Liu 2011; (12) Tamburri et al. 2012; (13) Eckmiller et al. 2011; (14) Diehl 
& Statler 2005; (15) Fukazawa et al. 2006; (16) Humphrey & Buote 2006; 
(17) Jetha et al. 2007; (18) Kim et al. 2007a; (19) Diehl & Statler 2007; (20) 
Rasmussen & Ponman 2007; (21) Jetha et al. 2008; (22) Jeltema et al. 2008; 
(23) Diehl & Statler 2008; (24) Wang et al. 2010; (25) Giacintucci et al. 2011; 
(26) Matsushita et al. 2012; (27) Crain et al. 2013; (28) Heida et al. 2013; (29) 
Gibson & Brandt 2012; (30) Ofek et al. 2013; (31) Immler et al. 2006; (32) 
Grier et al. 2011; (33) Smith et al. 2012. 


acis_f ind_afterglow. The pixel randomization normally 
used in the Chandra data pipeline was removed to prevent 
degradation of the spatial resolution. A O ! 5 pixel randomiza- 
tion is necessary for data with exposure times of <2 ks to 
compensate for aliasing effects, however, the exposure times of 
the observations in our sample are typically far in excess of this 
limit, therefore we omitted the randomization and recover the 
resolution to subtract robustly the point sources from the diffuse 
emission. 

For data taken in VFAINT mode, we applied the VFAINT 
cleaning algorithm, which uses a 5 x 5 pixel event island, rather 
than the 3x3 island in FAINT mode, for the rejection of 
cosmic rays. The VFAINT cleaning method has been shown to 
occasionally reject photons from bonafide X-ray point sources 
leading to underestimates of the associated fluxes, however, 
this does not impact our analysis of the diffuse emission. We 
then filtered the data on the standard Advanced Satellite for 
Cosmology and Astrophysics ( ASCA ) grades and selected only 
grades 0, 2, 3, 4, and 6 to produce the final Level 2 events file 
for analysis. 

Before examining the diffuse emission, we first excised the 
point sources from the observations. We used the Mexican hat 
wavelet detection routine wavdetect (Freeman et al. 2002) to 
search for point sources in the data. Multiple Chandra ObsIDs 
were merged using merge_obs prior to running the detection 
algorithm to facilitate the detection of faint point sources. The 
merged images were created after correcting the aspect solution 
to compensate for small offsets in the World Coordinate System 
between the multiple observations. 

To detect robustly all of the point sources in the data, we 
followed a prescription similar to that presented in Tiillmann 


Table 4 

Extraction Region Parameters 


Group 

J2000 Coordinates 
a 8 

Shape 

Radius 

HCG 30 

4 h 36 m 25 ! ;9 

-2°50' 14 5" 

Circular 

4'4 

HCG 37 

9 h 13 m 3ffi2 

29°59' 24 3" 

Circular 

3'3 

HCG 40 

9 h 38 m 55 : ;2 

-4°51' 2 1" 

Circular 

3.'8 

HCG 51 

H h 22 m 21?8 

24° 17' 39 5" 

Circular 

3'.5 

HCG 68 

13 h 53 m 36*7 

40° 18' 52 6" 

Circular 

6'3 

HCG 79 

15 h 59 m l 1?5 

20° 45' 26 2" 

Circular 

4.T 

HCG 97 

23 h 47 m 25*6 

-2° 19' 5 6" 

Elliptical 

3'.1 x 3'0 a 

HCG 100 

o h i m 2ao 

13°7' 2 8" 

Circular 

3'5 

RSCG 17 

l h 56 m 21 : ;6 

5°37' 53 6" 

Circular 

3'1 

RSCG 31 

9 h 17 m 23*4 

41°57' 17 7" 

Circular 

4'.3 


Note. a The position angle is 0° . 


et al. (2011). Specifically, we divided the data into the energy 
ranges 0.5-2 (“soft”), 2-8 (“hard”), and 0.5-8 keV (“full”) 
with block 1, 2, 4, and 8 pixel spatial binning in each energy 
range. This produced 12 images on which to run wavdetect. 
We generated a point-spread function (PSF) model for each 
position on the CCDs of interest using the mkpsf routine in 
CIAO with an encircled energy fraction of 95% at the midpoint 
of each energy range. The source significance threshold was 
dynamically set such that there was approximately one false 
source detected per wavelet scale in each image. Specifically, 
we used the falsesrc parameter in wavdetect to allow the 
source significance threshold for each pixel to vary. However, 
we note that individual, unbinned pixels cannot be used for 
source detection as wavdetect suppresses fluctuations on 
scales smaller than the PSF. We chose wavelet scales for source 
detection of 2 n / 2 for integers n such that 0 ^ n ^ 3. Our goal 
for the point source detection was to be extremely conservative 
and reject all potential point sources because, if present in 
the extracted spectra of the diffuse emission, they introduce 
strong biases in the results. We then matched the resulting point 
source catalogs using an angular separation tolerance of 075, 
first selecting the smallest spatial binning scale in which the 
source was detected in an energy band, ensuring the best centroid 
position for the source. These three catalogs were then matched 
with the same tolerance across the energy bands choosing the 
sources with the smallest PSF. This creates one point source 
catalog per CG. 

Extraction regions for the CGs were selected to include 
all of the member galaxies as well as any obvious diffuse 
X-ray emission. As discussed in Desjardins et al. (2013), the 
emission in most CGs is clearly not virialized, and therefore an 
extraction region defined according to the virial radius (e.g., 
rsoo) is not physically meaningful in these systems. Larger 
extraction regions simply result in additional noise and lead 
to larger uncertainties in our subsequent spectral-model fitting. 
Our extraction regions are defined for each group in Table 4. 

As in Desjardins et al. (2013), we used the ACIS stowed back- 
ground data for determination of the instrumental background. 
We note that in groups with low signal to noise, excess residual 
line emission was observed in the background- subtracted spec- 
tra at ~ 1.8 keV. A strong line at this energy is observed in the 
stowed background data, and we attribute the excess emission 
in the science spectra to under- subtraction of an instrumental 
feature, e.g., the Si K line at 1.845 keV or the iridium edge in 
the 1. 8-2.1 keV range. 
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Table 5 

Best-fitting X-Ray Model Parameters 



tbabs 

ztbabs 


MEKAL 




Group 

Galactic H i 
(10 20 cm -2 ) 

Redshifted H i 
(10 20 cm“ 2 ) 

kT 

(keV) 

z/z 0 

A a 

(10 -4 cm -5 ) 

Lx b 

(ergs -1 ) 

X 2 /dof 

HCG 30 
HCG 37 
HCG 40 
HCG 51 
HCG 68 
HCG 79 
HCG 97 
HCG 100 
RSCG 17 
RSCG 31 

4.72 

1.87 

3.60 
1.13 
0.96 
3.86 

3.61 
4.51 
4.33 
1.15 

io.i(£V£ 

18.57 ±3.14 

4 ga+1.31 
— 1.25 

0.6 

1 O2 +014 
0.6 

1.36 ±0.04 

0 57+0-07 

-0.08 

0.6 

0.85 ± 0.03 
0.6 

1 i c+0.04 
^ 1 — 0.03 

0.6 

0.5 

L21 ->\ 

0.5 

0.29 ± 0.04 
14.58+>\ 
0.5 

0.19 ±0.02 
0.5 
0.11 
0.5 

<0.47 
0 75+0-80 
7 —0.75 

<0.30 

29.067 1 9 8 4 3 

2-36f 05 24 

<0.28 

28.72^2.58 

<0.47 

21-78-Z23 

<0.24 

<40.59 
41.36f° 2 3 9 1 1 
<40.76 
42.73 ± 0.03 
41.13 ±0.06 
<40.54 
42.45 ± 0.04 
<40.66 
42.22 ± 0.02 
<39.70 

98.28/72 

100.97/77 

52.32/32 

311.61/145 

207.47/138 


Notes. Data for additional CGs can be found in Desjardins et al. (2013). 
a The normalization of the MEKAL model. 
b The X-ray luminosity over the range 0.01-100 keV. 

Events were extracted using the AC IS Extract (AE) software 
package 8 (Broos et al. 2010, 2012). The point source catalogs 
for each group were input into AE for PSF modeling using MARX 
version 4.4. The point sources were then excised from the events 
files prior to extracting the CG spectra. Specifically, we used AE 
to create a circular mask for each point source that enclosed 
99% of the PSF, and then multiplied the mask radius by a factor 
of 1.1 to ensure no contamination of the diffuse emission by 
the wings of the PSF. In addition, the point source masks were 
also applied to the stowed background data for the extraction 
of the background spectra. We used AE to generate response 
matrix files and the Cl AO tool mkwarf to create weighted 
ancillary response files using the weight map extension of the 
spectral files. 

Extracted spectra were then fit in XSPEC version 12.7.1 using 
a combination of foreground absorption and a thermal plasma. 
Prior to fitting, we binned the spectra using the HEASoft tool 
grppha such that each bin had a minimum of 20 counts; this 
ensures that Gaussian statistics (i.e., x 2 fitting) may be used. 
The multiplicative Tuebingen-Boulder interstellar medium ab- 
sorption model (tbabs) was used to account for photoelectric 
absorption along the line of sight. For this purpose, we used the 
relative abundances from Lodders (2003). The Galactic hydro- 
gen column density was fixed to the value from the weighted 
average of the Kalberla et al. (2005) H i maps using the HEASoft 
nH tool. The thermal plasma was modeled using the MEKAL 
plasma model (Mewe et al. 1985, 1986; Kaastra 1992, 1993; 
Liedahl et al. 1995; Kaastra & Liedahl 1995) with the ioniza- 
tion balance from Arnaud & Rothenflug (1985) and Amaud & 
Raymond (1992). The low energy resolution non-grating spectra 
are insufficient to determine plasma densities, therefore we fixed 
the density in the model to a reasonable value of n = 1 cm -3 . We 
calculated the X-ray luminosities over the range 0.01-100 keV 
using a “dummy” response created by the XSPEC command 
dummyrsp. The results of the best-fitting spectral models are 
reported in Table 5. 

In HCGs 51 and 97, and RSCG 17, the model overestimated 
the X-ray emission below ~0.7 keV. We used HCG 97 to 
test three different additional model components: (1) a second 
MEKAL plasma; (2) a simple power law; and (3) additional 


8 The AC IS Extract software package and User’s Guide are available at 

http://www.astro.psu.edu/xray/acis/acis_analysis.html. 


absorption at the redshift of the groups (ztbabs in XSPEC). The 
second plasma component of case (1) failed to fit the observed 
flux at low energies; however, the power law and the additional 
absorption produced nearly equal values of x 2 = 2.0 and 2.1, 
respectively. In both cases (2) and (3), the temperatures were 
identical within the errors. In case (2), the power law had a 
hard photon index of 1.4, which led to a higher luminosity 
compared to the additional absorption model (log 10 [Lx, p i] — 
log 10 [Lx, ztbabs] = 0.41) due to increased flux at higher energies. 
We chose to use the additional absorption component and 
find redshifted Hi column densities in HCGs 51 and 97, and 
RSCG 17 of 1.01 x 10 21 , 1.86 x 10 21 , and 4.93 x 10 20 cm" 2 , 
respectively, corresponding to H i masses of ~10 5 M©. We note 
that HCGs 5 1 and 97 have very extended X-ray emission (see 
Section 4.1), and if the absorption interpretation is correct, we 
may be detecting low-surface brightness cool gas on the near 
sides of these systems. Indeed, this would be consistent with 
the H i upper-limit of HCG 97 with Very Large Array L-band 
imaging (S. Borthakur 2013, private communication). The X-ray 
emission in RSCG 17 subtends a much smaller angle compared 
to HCG 97, but may still be explained with the low-surface 
brightness interpretation. Comparatively, the inclusion of a hard 
power-law component does not have a physical motivation, but 
cannot be summarily ruled out because of the low signal at 
E> 3 keV. 

If the number of X-ray photons associated with the CG was 
less than 3a above the instrumental background, we classified 
such a source as a non-detection and used a plasma temperature 
of T = 0.6 keV to set an upper-limit on the X-ray luminosity, 
identical to the method presented in Desjardins et al. (2013). 
We use the definition of a from Desjardins et al. (2013) such 
that a = [SB + ( A s t s / A b tb)B ] l/ 2 , where SB is the total counts 
in the source before background subtraction, B is the number of 
counts in the background, A is the area of the extraction region, 
t is the integration time, and the subscripts a and b represent the 
science and background observations, respectively. 

Note that HCG 79, also known as Seyfert’s Sextet, lies at 
Galactic coordinates i = 35?0 and b = 46? 9 and is coincident 
with a portion of the North Polar Spur (NPS). The NPS is a 
region of bright, soft X-ray emission associated with expanding 
supernovae remnants (e.g., Cruddace et al. 1976; Borken & I wan 
1977; Iwan 1980; Miller et al. 2008), therefore making it more 
difficult to detect emission from HCG 79. In addition, an X-ray 
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bright background group or cluster of galaxies with z ~ 0.3 is 
located (X6 to the northwest of HCG 79 (Palma et al. 2002; 
Tamburri et al. 2012). Rather than spatially model and exclude 
the emission from this background source, we included it in 
the spectral extraction and then modeled it with an additional 
plasma component. 

3.2. Optical Data 

To compare the diffuse X-ray and the relative, optical bright- 
nesses of the two brightest group galaxies in each CG (see 
Section 4.3), we used observations of HCGs 16, 19, 26, 33, 40, 
42, 48, and 62, and RSCG 15, obtained on 2011 January 18 at 
the APO 3.5 m telescope using the Seaver Prototype Imaging 
camera (SPIcam) instrument and the SDSS r' filter. SPIcam is a 
2048 x 2048 pixel CCD with a scale of O'.' 14 pixel -1 ; however, 
because the APO site is seeing-limited, we used 2x2 pixel bin- 
ning to facilitate faster readout time without sacrificing spatial 
resolution. For the remaining CGs in the Walker et al. (2012) Ex- 
panded Sample, we used r'-band images taken from the SDSS 
DR9 database. We include all of the CGs in the Expanded Sam- 
ple rather than only the CGs observed with Chandra to have a 
statistically large enough sample to compare against the X-ray 
groups. Note that between the APO and SDSS observations, 
we have optical coverage of the entire Expanded Sample except 
HCGs 90 and 91, and RSCG 4. This left us with stellar mass 
measurements for galaxies in 47 CGs. The optical photometry 
of CGs using SDSS data is further explored in Walker et al. 
(2013). 

We reduced the APO r'-band data using PyRAF version 2.0 
and I RAF 9 version 2.14 to perform serial overscan subtraction 
and to create master two-dimensional bias and dark frames, as 
well as master flat images in both filters. We used Source Extrac- 
tor (Bertin & Arnouts 1996) version 2.8.6 for all extended source 
photometry with an aperture set to twice the Petrosian radius to 
ensure uniformity across the sample. No absolute photometric 
calibration was performed as we were only interested in the 
differential photometry of the two brightest group galaxies. The 
small projected separations of galaxies in CGs necessitated that 
we be particularly careful with object blending. In cases where 
multiple sources overlap, Source Extractor uses a de-blending 
algorithm to separate the pixels associated with each object. We 
found that the default de-blending parameters were sufficient 
for E/ SO galaxies; however, we needed to adjust the settings on 
an individual basis for inclined, star-forming galaxies to ensure 
that the entire galaxy was classified as one source rather than a 
collection of blended sources. 

Combined with the SDSS DR9 images, we have optical 
photometry in the r' band for 41 CGs, all but one of which are 
in the Walker et al. (2012) Expanded Sample, while HCG 51 is 
solely in our X-ray sample. Note that HCG 30, which is in our 
X-ray sample, but not the Walker et al. (2012) Expanded Sample, 
does not have APO or SDSS data. We use the Hickson (1982) 
ordering of CG galaxies according to their optical brightness 
to select the first and second rank galaxies with respect to the 
POSS £-band luminosity. For the RSCGs, we use the ordering 
presented by Walker et al. (2013), which uses the same ordering 
system albeit in a marginally different bandpass compared to 
the HCGs. 


9 IRAF is distributed by the National Optical Astronomy Observatory, which 

is operated by the Association of Universities for Research in Astronomy 
(AURA) under cooperative agreement with the National Science Foundation. 


3.3. Stellar Mass Determination 

The standard method of stellar mass determination, i.e., use 
of the A^-band luminosity, assumes a universal mass-to-light 
ratio independent of the galaxy morphology and, therefore, star 
formation history. To determine more robustly the stellar masses 
of the galaxies in our sample, we used the library of galaxy 
templates generated by the GRASIL 10 code (Silva et al. 1998, 
2003; Silva 1999, 2009; Granato et al. 2000; Bressan et al. 2002; 
Panuzzo et al. 2003; Vega et al. 2005) to fit the galaxy spectral 
energy distributions (SEDs). Two Micron All Sky Survey 
(2MASS) JHK S and Spitzer Infrared Array Camera (IRAC) 
3. 6-8.0 /im fluxes were taken from Walker et al. (2012). In the 
case of saturation in one of the four IRAC bands (six galaxies), 
or if IRAC data were missing (nine galaxies), we calculated 
the stellar mass using only the A^-band M/L relation from 
Bell et al. (2003), i.e., M 0 /L y 0 = 0.95 zb 0.03. For galaxies 
missing 2MASS data from Walker et al. (2012), we used the 
2MASS photometry from the WISE database. For all sources, we 
converted the fluxes to luminosities using distances determined 
from the 3 K cosmic microwave background (CMB) dipole- 
corrected velocities. All of the elliptical and spiral templates, 
as well as the M82 starburst galaxy template, were fit to the 
data without knowing a priori the galaxy morphology to ensure 
unbiased results. We note that the best- fitting templates do agree 
well with the observed galaxy morphologies in a general sense, 
i.e., spiral galaxies are best modeled using spiral templates and 
likewise for elliptical galaxies, though the finer divisions within 
these classes (e.g., Sa, Sb, and Sc) are sometimes not accurately 
determined from the SED fitting. 

In the near-infrared, the SED of a galaxy scales with stellar 
mass, therefore the normalization of the best-fit galaxy template 
to the observed luminosity coupled with the stellar mass of the 
model yields the stellar mass of the galaxy. To properly fit the 
templates to the galaxy photometry, we first shifted the GRASIL 
SED templates to the observed frame of the source and then 
convolved them with the 2MASS and Spitzer filter response 
curves. We used the normalization to the 3.6 pm luminosity 
as an initial guess of the stellar mass normalization before 
X 2 minimization. We note that mid-IR active galaxies contain 
strong polycyclic aromatic hydrocarbon emission features in 
the 5.8 and 8.0 pm IRAC bands; however, these did not greatly 
affect the quality of the SED fitting. We estimated the errors 
on the SED-fitted masses by varying the normalization of the 
template until A/ 2 = 2.71, i.e., the 90% confidence interval. 
The total stellar masses for the CGs in Desjardins et al. (2013) 
and this paper are listed in Table 6. The masses of the individual 
galaxies in the Expanded Sample are shown in Figure 1 where 
we compare the SED fitted masses against the stellar masses 
derived from the X s -band M/L relationship from Bell et al. 
(2003). We find that the stellar masses from SED fitting match 
well with some scatter compared to those from the A^-band 
method, though there is a small deviation at low masses. 

4. RESULTS AND DISCUSSION 

We detect diffuse X-ray emission in 50% of the CGs in our 
new sample observed with Chandra. Combined with the CGs 
from Desjardins et al. (2013), this yields 19 groups with 12 
detections and an overall detection rate of 63%. We caution 
the reader that our detection rate should not be used to draw 
conclusions about the statistical distribution of diffuse X-ray 


10 http://adlibitum.oats.inaf.it/silva/grasil/grasil.html 
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Table 6 

Total Group Stellar Masses Using Core Galaxies 


Group 

Stellar Mass 
(10 9 Mo) 

Ngal 

HCG 7 

334.39 ± 0.34 

4 

HCG 16 

462.24 ± 1.92 

4 

HCG 22 

104.58 ± 0.23 

3 

HCG 30 

0.88 ± 0.00 

4 

HCG 31 

32.37 ± 0.26 

3 

HCG 37 

460.67 ± 0.28 

5 

HCG 40 

398.47 ± 0.40 

5 

HCG 42 

458.91 ±0.36 

4 

HCG 51 

49.68 ± 0.00 

5 

HCG 59 

36.05 ± 0.00 

4 

HCG 62 

229.34 ± 0.37 

4 

HCG 68 

427.37 ± 0.56 

5 

HCG 79 

79.56 ± 0.01 

4 

HCG 90 

308.05 ± 1.61 

4 

HCG 92 

523.66 ±2.13 

4 

HCG 97 

328.03 ± 0.44 

5 

HCG 100 

127.11 ±0.00 

4 

RSCG 17 

288.40 ± 2.30 

3 

RSCG 31 

72.96 ± 0.07 

3 


luminosities of CGs (e.g., Ponman et al. 1996) as many of 
the targets that make up our sample were observed on the 
assumption that the groups would be X-ray bright. 

4.1. X-Ray Morphology 

Desjardins et al. (2013) find that, in contrast to galaxy 
clusters, the diffuse X-ray emission in CGs is often linked to the 
individual galaxies rather than the group itself. Therefore, we 
construct contour maps of the X-ray emission to examine the 
distribution in this extended X-ray sample. To make the contour 
maps, we first excised the point sources and interpolated over 
them using the CIA0 task dmf ilth. Note that these interpolated 
images were only used in the creation of the contour maps 
and not in the spectral analysis of the diffuse X-ray emission. 
The resulting image of only diffuse emission was divided by 
the monoenergy exposure map (optimized at a photon energy 
of 1 keV) to create flux images in units of photons s -1 cm -2 . 
Finally, we smoothed the flux images by convolving them with 
a Gaussian kernel. 

Figure 2 shows the diffuse X-ray contour maps for the X-ray 
detected CGs (see Desjardins et al. 2013 for the contour maps 
of the CGs in that paper). We once again find a mixture of 
galaxy- and group-linked emission, with HCGs 5 1 and 97 and 
RSCG 17 having the most extended X-ray halos. In all cases, 
detected X-ray emission is centered on the optically brightest 
group galaxy, and is not as localized as in HCGs 16 and 31 
(Desjardins et al. 2013). This suggests that for the X-ray detected 
CGs with galaxy-linked diffuse emission in this paper that were 
not previously presented in Desjardins et al. (2013), the sources 
are hot gas halos around the brightest group galaxies and not 
star formation within the galaxies. This is further supported by 
the fact that all of the X-ray detected CGs in this new sample of 
10 groups have E/S0 brightest group galaxies. 

In HCGs 68 and 97, the X-ray morphologies are indicative 
of recent or ongoing galaxy-galaxy interactions. In Figure 2, 
one can see a hot-gas bridge connecting HCG 68 A and B. 
Interestingly, the morphology of the X-ray gas in HCG 97 is 
elongated in three directions away from the brightest group 
galaxy toward galaxies 97D and E, and is most pronounced 



M/L Mass log 10 (M*) [M 0 ] 

Figure 1 . Comparison of the stellar masses determined by mid-IR SED fitting 
against the masses using the K- band M/L relationship from Bell et al. (2003). 
Galaxies in this X-ray study are plotted with filled symbols while the remaining 
galaxies from the Expanded Sample defined by Walker et al. (2012) are plotted 
as open symbols for comparison. The dashed line in the upper panel shows the 
one-to-one relation for reference, while the solid line is the orthogonal distance 
regression fit to the data. The bottom panel shows the residuals around the fit 
with a dashed line for reference. Most of the observed scatter is likely due to 
scatter in the M/L relation, as well as uncertainties in the SED fitting. 

(A color version of this figure is available in the online journal.) 


toward the southeast. This last direction may indicate a past 
tidal encounter with one of the other galaxies in the group. 


4.2. X-Ray Scaling Relations 

We wish to investigate how the CG X-ray properties compare 
with the X-ray cluster scaling relations between the bolometric 
X-ray luminosity ( L x ), X-ray temperature, and group velocity 
dispersion. This allows us to examine the physical nature of 
the hot plasma in groups independent of the gas morphology. 
These power-law scaling relations are expected from the virial 
theorem and the self- similar model of galaxy cluster and ICM 
formation (Kaiser 1986). 

After inspecting the hot gas morphologies and surface bright- 
ness profiles of the CGs, we surmise that the six most luminous 
CGs may have additional X-ray emission out to larger radii (see 
HCG 62 in Desjardins et al. 2013). To compensate for this, we 
used a (hydrostatic, isothermal sphere) model to quantify the 
X-ray luminosity correction out to r$ oo- The /3 model of the 
surface brightness profile is given as 


S(R) = So 



2-1 — 3 / 3 + 0. 5 


( 1 ) 


where So is the peak surface brightness, R is the distance from 
the center, and r c is the core radius (for which we have once 
again assumed the median two-galaxy separation). The /3 term 
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Right Ascension [deg] 


Figure 2. X-ray contour maps of the CGs in the 0.5-2 keV band. The optical images come from the SDSS DR9 and are in the r' filter. In HCGs 37, 68, 97, and RSCG 
17, the X-ray emission is marked by green contours at levels of 2.5 x 10 -8 , 5 x 10 -8 , 1 x 10 -7 , and 2.5 x 10 -7 photons s -1 cm -2 . In HCG 51, the green X-ray 
contours correspond to 1 x 10 -7 , 1.5 x 10 -7 , 2 x 10 -7 , and 2.5 x 10 -7 photons s -1 cm -2 . The red dashed region indicates the outer extraction boundary for each 
CG. X-ray point sources are also excised, but are not labeled in these images. We label CG members with accordant redshifts using the notation from Hickson (1982) 
and Walker et al. (2013) for reference. 

(A color version of this figure is available in the online journal.) 


is defined as 


p = 


/jim p (7 2 

kT ’ 


( 2 ) 


where /x is the mean molecular weight (fixed at solar), m p is the 
mass of the proton, a is the line-of-sight velocity dispersion, k 


is the Boltzmann constant, and T is the gas temperature. This 
resulted in correction factors of the order of unity in most cases 
except HCG 62, which required a factor of 12.8 to compensate 
for the small extraction region used in Desjardins et al. (2013). 
Note that in Desjardins et al. (2013), we found a correction 
factor of 3.1 was necessary to scale the observed Chandra flux 
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X-ray Temperature log 10 (T) [keV] Velocity Dispersion log 10 (cr) [km s" 1 ] X-ray Temperature log 10 (T) [keV] 


Figure 3. X-ray scaling relations for the cluster samples of Wu et al. (1999) and Zhang et al. (2011), corrected to our cosmology, and the CGs: Lx~T (left); L x -cr 
(center); and cr-T (right). In all three panels, the gray crosses represent the galaxy clusters, while the filled and open circles are the X-ray detected and non-detected 
CGs, respectively, in the combined sample of this study and Desjardins et al. (2013). In the Lx~T and Lx~cr diagrams, the dashed line shows the orthogonal distance 
regression fit to the clusters. Such a fit is not possible for the cr-T relation because of the large scatter in the cluster data and the paucity of clusters at low velocity 
dispersions. For the non-detected CGs, we use a reasonable temperature of 0.6 keV for determining upper-limits in Lx. 


to that observed by ROSAT in the much larger aperture used by 
Mulchaey & Zabludoff (1998). 

The left panel of Figure 3 shows the CGs and the galaxy 
cluster sample from the literature in the L x -T plane. There 
exists a population of high-temperature groups (T ~ 1 keV; 
HCGs 51, 62, 97, and RSCG 17) that agree well with the 
cluster L x -T relation, and two lower temperature CGs (HCGs 
37 and 42) that are also in agreement. The agreement between 
the aforementioned hot, X-ray luminous CGs and the clusters 
occurs within the scatter of the cluster data. In Desjardins et al. 
(2013), the authors find groups that agree with the cluster scaling 
relations are those in which the emission is linked primarily to 
the IGM rather than to the individual galaxies (e.g., HCG 62); 
however, HCG 37 represents a CG where the ionized gas is 
clearly associated with the brightest group galaxy, and therefore 
may be probing the cool- temperature, low-luminosity portion 
of the X-ray scaling relations. If this is true, it may indicate that 
HCG 42, which has a similar temperature and X-ray luminosity, 
is not an example of group-linked emission, but a CG in which 
the X-ray halo is only observed around the brightest group 
galaxy. 

The CGs and clusters in the L x -o plane are shown in the 
center panel of Figure 3. Consistent with the L x -T relation, 
there exists a population of high velocity dispersion CGs which 
seem to agree well, albeit with more scatter, with the cluster 
sample from the literature, and in fact these are the same CGs as 
those that matched the L x -T relation. These six CGs, HCGs 37, 
42, 51, 62, 97, and RSCG 17, may represent the most cluster- 
like CGs in our sample with respect to their hot gas properties. 
We again note that the consistency between the high-dispersion, 
X-ray luminous CGs and the clusters occurs at the level of 
the scatter in the cluster data, while the CGs themselves all 
exist systematically below the L x -o relation fitted to the galaxy 
clusters. 

For completeness, we also examine the cr-T relationship for 
CGs and clusters in the right panel of Figure 3, however, we 
find that the clusters exhibit a large scatter and, due to under- 
sampling at lower velocity dispersions/cooler temperatures, the 
relationship is too poorly constrained for comparison with the 
CGs. We do note that there is a rapid drop in velocity dispersion 
at cool X-ray temperatures in the CGs, further indicating that 
these are not relaxed systems. 


4.3. Relating X-Ray Emission to the Baryonic Mass 

Understanding the relative importance of the hot, warm, and 
cool gas phases is critical to understanding the evolution of 
systems of galaxies. At a more fundamental level, the distri- 
bution and phase of the baryons in galaxy groups dictate the 
future evolutionary path of the system, while also giving in- 
sight into the group history. For example, a system lacking 
cool/cold gas while having a hot X-ray halo (e.g., HCG 62; 
Desjardins et al. 2013) is unlikely to convert much more gas 
into stellar mass, and therefore future galaxy evolution will be 
primarily dynamical. Conversely, galaxies in systems such as 
HCG 16 that are Hi-rich will evolve both dynamically and in 
terms of their stellar populations. While the reservoir of cool 
and cold gas in systems of galaxies is critical to these two ex- 
amples, it is likely that the bulk of the baryonic mass will end 
up in either stars or an X-ray emitting halo, 11 therefore exam- 
ining the relationship between stellar mass, cool gas, and group 
X-ray emission gives a more complete picture of galaxy evolu- 
tion in groups. We hypothesize that the mass of these systems is 
critical to how they will evolve, i.e., it will determine the ability 
of the group potential to heat gas to X-ray temperatures and 
consequently lower the baryon fraction in stars. An increase in 
the amount of hot gas between the galaxies will further affect 
galaxy evolution within the groups through, e.g., ram-pressure 
stripping. 

If we disregard the dynamical masses due to their large 
intrinsic uncertainties (see Section 4.4), then we can compare 
the relative masses of the groups as the sum of their directly 
observable components. For example, HCG 22 has no detected 
diffuse X-ray emission (log 10 [Lx] < 39.84 erg s _1 ; Desjardins 
et al. 2013), therefore we can assume that relatively few baryons 
are found in the hot IGM. Furthermore, the molecular gas mass 
of non-star-bursting galaxies is typically negligible compared 
to the combined Hi and stellar mass (see, e.g., Young & 
Scoville 1991), thus the total baryonic mass of the group may be 
approximated as the sum of the stellar and H i masses (similar 
to Connelly et al. 2012). If we assume that the fraction of 
group mass in baryons, as in more massive clusters, is constant 

1 1 This neglects the contribution of other negligible components of baryonic 
mass in systems of galaxies (e.g., dust), as well as baryons that are expelled 
from groups. 
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Figure 4. CG diffuse X-ray luminosity as a function of the total stellar and H i 
mass. Assuming that the baryon fraction of groups is approximately constant, 
as has been seen in clusters, the total stellar and Hi masses then trace the 
total masses of low-mass systems, while they provide lower-limits on the total 
masses of X-ray luminous groups with larger hot gas masses. The lower-limits 
in the total mass arise from the exclusion of additional group members at large 
distances from the compact core region, but which are still within r$ oo- 


(e.g., Gonzalez et al. 2007; Andreon 2010), then HCG 22 is at 
least a factor of four less massive than the X-ray luminous HCG 
42. This example initially supports our hypothesis regarding the 
existence of an X-ray halo and its dependence on the group 
mass, therefore we choose to expand our test to a larger sample. 

We plot the diffuse X-ray luminosity against the total stellar 
and H i masses in Figure 4. While the most X-ray luminous and 
cluster-like CGs likely have a substantial fraction of baryons in 
the hot phase, it is not possible to measure the hot gas mass 
of baryons without making potentially poor assumptions about 
the diffuse emission, e.g., spherical symmetry and the value of 
dT /dr , therefore we do not include the hot gas mass in our 
total mass estimates. Again assuming that the baryon fraction in 
groups is approximately constant, then the total stellar and H i 
mass may be used as a proxy for the total group mass in X-ray 
faint systems, and as a lower- limit in X-ray luminous CGs. From 
Figure 4, we find that nearly all of the X-ray luminous CGs with 
intragroup X-ray hot gas have higher total stellar and H i masses 
than the CGs that were not detected by Chandra and those 
with galaxy-linked X-ray emission. The approximately vertical 
distribution of CGs at M* + M H i =11.6 M 0 is likely due to an 
increasingly large mass of baryons in the X-ray phase in hotter, 
more X-ray luminous groups. Indeed, these groups have line-of- 
sight velocity dispersions that imply relatively high dynamical 
masses, and therefore a substantial fraction of baryons may be 
in the hot phase. The three low-mass CGs with detectable X-ray 
emission are HCGs 31, 51, and 59, two of which have X-ray 
emission associated with vigorous star formation (Desjardins 
et al. 2013). Excepting HCG 51, these results support our 
hypothesis that low-mass groups do not heat their IGM to high 
temperatures. We caution the reader that the hot gas in low-mass 
systems may have very high entropy, and therefore low density. 
From the proportionality L x oc p 2 , we expect that such gas 
would be very difficult to detect. Therefore, HCG 51 may in 
fact have a lower total mass than the others, or it may have a 
higher fraction of its mass in the hot gas phase (i.e., the gas may 
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be low entropy as discussed above). Finally, we note that there 
are two groups, HCGs 7 and 40, with high total stellar and Hi 
masses but no detectable X-ray emission. 

We note that previous work has examined loose groups of 
galaxies without detectable X-ray emission, and the consensus 
is that the gas in these groups is too cool to produce significant 
X-ray luminosity (Mulchaey et al. 1996; Rasmussen et al. 2006). 
The lack of hot gas has been explained in two different ways: 
Mulchaey et al. (1996) argue that some groups are too low- 
mass to effectively heat their gas to X-ray temperatures; while 
Rasmussen et al. (2006) hypothesize that X-ray underluminous 
groups are dynamically young (i.e., in the process of collapsing) 
and have not had sufficient time to virialize their IGM. These 
ideas are not mutually exclusive, and indeed both may affect 
the formation of group X-ray halos. Further, we hypothesize 
that the absence of an observable hot IGM in low-mass groups 
may be caused by the hot gas having very low density, in which 
case it would be have very low surface brightness and be nearly 
impossible to detect, or that the gas is never heated to very high 
temperatures and instead cools efficiently. Note that our two 
hypotheses are independent from one another as the gas density 
must be high for efficient cooling to occur. Unfortunately, our 
data are insufficient to test the low-density gas and efficient 
cooling scenarios. 

Regarding the explanations offered by Mulchaey et al. (1996) 
and Rasmussen et al. (2006), we present the two examples of 
HCGs 31 and 7. In HCG 31, we find an unusually high baryon 
fraction of ft ~ 0.36 (assuming negligible mass in the form of 
hot gas) due to the high H i mass, while the ratio of stellar mass 
to dynamical mass is relatively large, and therefore may indicate 
that the dynamical mass is underestimated. We remind the reader 
that the dynamical mass is very uncertain due to the small 
number of galaxies available with which to measure the velocity 
dispersion, and the magnitude of this uncertainty is unclear. 
Assuming that most of the baryonic mass in HCG 31 is in the 
stellar and H i components, then a reasonable baryon fraction 
implies it is a relatively low-mass system. Thus, we expect HCG 
3 1 to have a very cool virial temperature. Conversely, HCG 7 is 
a relatively massive group (M dyn = 10 121 M©; Desjardins et al. 
2013), and therefore we expect a hotter virial temperature based 
solely upon the total mass. Again, we regard the dynamical mass 
with caution, though we do note that its stellar mass implies that 
it is at least an order of magnitude more massive than HCG 3 1 . 
Thus, HCG 7 seems X-ray underluminous for its mass, though it 
is not dynamically young. Indeed, Konstantopoulos et al. (2010) 
find that the galaxies in HCG 7 show evidence for long-term, 
enhanced evolution in the group environment without direct, 
strong interactions. It is unclear how HCG 7 may form an X-ray 
halo at some later evolutionary stage, if it will at all. 

Jones et al. (2003) observationally define fossil groups to 
identify groups near the end of their evolution in which the 
groups are dominated by a single galaxy and have X-ray 
emission in excess of that associated with normal galaxies. The 
authors use as their criteria a diffuse X-ray luminosity above a 
specific threshold ( L x > 10 42 h ^ erg s -1 ) and a difference in 
R-band magnitude between the two brightest group members 
Am 12 > 2 mag; however, the Am 12 criterion assumes that the 
group luminosity function is a well- sampled Schechter function 
(Schechter 1976), i.e., that the brightest group galaxy is much 
more luminous than L*, and therefore is located at the bright 
end of the exponential portion of the luminosity distribution. 
The Hickson (1982) CG selection criteria require that group 
members have Am ^ 3 mag with respect to the brightest group 
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Figure 5. Fraction of the total CG stellar mass contained within the most massive group member as a function of the SDSS r'-magnitude difference between the 
two brightest CG galaxies left) and the stellar mass ratio between the two most massive group galaxies (right). Groups that have an X-ray luminosity 

log 10 (Lx) ^ 41.4 erg s -1 (from the Jones et al. 2003 fossil group criteria corrected to our cosmology) are marked with diamonds. One of the magnitude differences 
in the left panel is an upper-limit because the APO images of HCG 42 A and 42B were taken at two difference airmasses. The r-band magnitude difference appears 
to be weakly correlated to the ratio of stellar mass contained within the most massive group galaxy compared to the rest of the group. Therefore, we suggest that the 
difference in stellar mass between the two most massive galaxies is a better indication of the evolutionary state of the group. The magnitude difference will be related 
to the stellar mass ratio assuming a single M/L ratio, but this assumption is not justified given the diversity of galaxies found in the most massive CG galaxies. 
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galaxy, making such a luminosity distribution unlikely, and 
therefore the Jones et al. (2003) criteria would have included 
no HCGs as fossil groups. We note that Dariush et al. (2010) 
use a magnitude difference of Am 14 > 2 mag, thus identifying 
more robustly low-mass fossil systems and potentially several 
CGs; however, the optical magnitude selection criterion of the 
HCG catalog makes it likely that this definition would still miss 
several fossil groups. Indeed, only HCG 42 comes close to 
satisfying the Jones et al. (2003) fossil group selection criteria. 
We thus seek to form an alternate physically motivated definition 
of fossil groups in order to classify the CGs in our X-ray sample. 

In Figure 5, we compare the optical selection criterion from 
Jones et al. (2003) against the distribution of stellar mass within 
the CGs in the Walker et al. (2012) Expanded Sample. 12 We 
find that the optical magnitude difference Am 12 between the 
brightest galaxies does not trace well the concentration of stellar 
mass in the CGs, and instead suggest that the fraction of the 
group stellar mass contained within the first-ranked galaxy (i.e., 
the group member with the largest stellar mass) can be better 
determined using the difference in stellar mass between the two 
most massive group galaxies. Note that we strictly use only the 
two most massive galaxies in the stellar mass criterion, and only 
the two brightest galaxies in the SDSS r' filter for calculating 
Am 12. In groups with two dominant galaxies approximately 
equal in mass, the most massive and brightest galaxy may not 
be the same; though, this is only possible when both Am 12 and 
log 10 (Mi/M 2 ) are small. 

The existence of a lone, massive, early-type galaxy has three 
possible implications for the groups: that at least one major 
merger has resulted in the formation of an elliptical galaxy and 
that only minor mergers will occur in the future; a series of 
minor mergers has concentrated the bulk of the stellar mass into 
a cold-gas-poor lenticular galaxy; or there was only ever one 


12 We excluded RSCGs 67 and 68 from this analysis as they are comprised of 
galaxies in the core of the Coma Cluster, as well as RSCG 32, which is 
embedded within the A779 cluster. 


massive galaxy in the group and it is now cold-gas-poor. In all 
three cases, the bulk of the stellar mass exists in a “red and 
dead” galaxy that has reached the end of its evolution in most 
respects. We therefore use the stellar mass distribution in our 
classification of CGs as evolved fossil systems. Specifically, we 
require that the first-ranked galaxy contain ^60% of the group 
stellar mass; or, alternatively, the first-ranked galaxy be at least a 
factor of three more massive than any other group member. We 
find that groups evolve from the lower-left region of the right 
panel of Figure 5, in which the two most massive galaxies are 
of approximately equal mass, to the upper-right region, where 
groups are dominated by a single massive galaxy. The values 
used in our fossil group selection criteria (i.e., 60% of the group 
stellar mass or a factor of three difference in the stellar masses 
of the first- and second-ranked galaxies) reflect the distribution 
of CGs in the right panel of Figure 5. Rather than require 
both of these criteria to be true, we only require that one be 
satisfied to be considered a fossil group candidate due to the 
scatter in the CG distribution; in either case, the stellar mass is 
clearly concentrated in a single group member. To ensure that 
groups we classify as fossil groups are indeed highly evolved, 
we also impose a morphology criterion requiring that the first- 
ranked galaxy be a E/S0 galaxy. 

We further remove the X-ray criterion from Jones et al. (2003) 
because it assumes that the potential of the group is sufficient to 
virialize the gas to hot temperatures, which is not true in low- 
mass groups. Indeed, if the gas in low-mass groups is heated 
to T < 10 6 K, then it will cool quickly and may therefore be 
available again for star formation (Dalgamo & McCray 1972; 
Sutherland & Dopita 1993; Schure et al. 2009). We note that 
the X-ray luminosity requirement from Jones et al. (2003) was 
intended to ensure that only truly bound systems were classified 
as fossil groups. Regardless of whether or not the systems 
are bound, the location of the galaxies in mid-IR color-color 
space shows evidence for accelerated evolution not observed in 
other environments (Johnson et al. 2007; Walker et al. 2010, 
2012), therefore the galaxies are in close physical proximity 
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for long enough to affect their evolution in measurable ways. 
Furthermore, the isolation criterion imposed on the HCGs by 
Hickson (1982) makes it unlikely that there are other galaxies 
to which the CG members may be bound. 

Finally, 21 CGs in the Expanded Sample from Walker et al. 
(2012) meet our first two fossil group criteria concerning the 
distribution of stellar mass; however, only six of these CGs 
(HCGs 19, 22, 40, and 42, and RSCGs 44 and 86) have 
an E/SO first-ranked galaxy and are therefore fossil groups 
under our new definition. From the three fossil CGs with 
Chandra observations, only HCG 42 has a substantial X-ray halo 
according to the Jones et al. (2003) definition, which requires 
a clarification of previous findings that fossil groups agree 
with the galaxy cluster X-ray scaling relations. Our findings 
suggest that it is only the massive, X-ray luminous fossil groups 
that agree well with the cluster scaling relations, while low- 
mass fossil groups systematically fall below these relationships. 
Furthermore, there may be many fossil groups that have not been 
classified as such because they are not massive enough to host 
X-ray luminous halos. Thus, if the X-ray luminous fossil groups 
represent the origin of optically and X-ray bright, isolated field 
ellipticals (e.g., Mulchaey & Zabludoff 1999), then fossil groups 
with no X-ray emission may be the progenitors of low-mass 
isolated ellipticals. Indeed, Mulchaey & Jeltema (2010) find 
that low-mass field ellipticals are typically X-ray underluminous 
and suggest that such galaxies, in contrast to their massive 
counterparts, may not be able to retain hot gas halos primarily 
due to strong winds from supernovae and, as a secondary 
factor, active galactic nucleus. In contrast, O’ Sullivan & Ponman 
(2004) argue that while the galaxy mass may be important 
in retaining X-ray halos, they hypothesize that these galaxies 
must also be the result of recent mergers to have sufficiently 
strong supemovae winds that can drive away low-density hot 
gas; though the authors also find that the metallicities of field 
ellipticals do not support supernovae-driven wind models. 

4.4. X-Ray Emission and the Hi Reservoir 

In Desjardins et al. (2013), the authors compare the 
X-ray luminosity to the Hi mass normalized to the dynami- 
cal mass and find tentative evidence to support the hypothesis 
by Konstantopoulos et al. (2010) that the X-ray emission in CGs 
may be dependent upon the morphology of the Hi gas in the 
system. Systems in which the H i has been stripped out of the 
galaxy disks into the IGM have material to fuel an X-ray lumi- 
nous halo. We note that the H i stripped from the galaxy disks 
may not be the only source of fuel for the hot IGM, however, 
a full analysis of the origin of such gas is beyond the scope 
of this paper. The explanation regarding the morphology of the 
Hi now seems secondary to the mass of the group given our 
comparison of the X-ray luminosity and the total stellar and 
Hi mass in Section 4.3. For example, the mass of the group, 
and the individual masses of the group members, will cause 
differences in the distribution of gas in the intragroup space. 
However, we still compare the X-ray luminosity to the H i reser- 
voir to understand better how neutral gas is consumed in CGs. 
A concern when using the velocity dispersions, and thus the dy- 
namical masses (as in Johnson et al. 2007 and Desjardins et al. 
2013), is the inherent uncertainty that stems from the small 
population of galaxies in each CG. Therefore, to mitigate this 
uncertainty, we compared the diffuse X-ray luminosity to the 
ratio of the Hi mass to the total stellar and Hi mass (fb,m) in 
Figure 6. 



lo(^ b,Hl) 

Figure 6. Diffuse X-ray luminosity as a function of the Hi mass relative to 
the total stellar and Hi mass ( fb,m ). Filled circles represent CGs with a total 
stellar and H i mass of > 1 0 1 1 3 Mq , while empty diamonds are CGs with a total 
stellar and Hi mass of <10 113 Mq. The lower limits on the values of fb,H\ 
represent the groups with additional members that may cause underestimates in 
the stellar mass measurements. The relatively massive CGs (i.e., those shown 
as large cirlces) are more luminous in X-rays and extend to smaller values of 
fb, Hi, and therefore more gas in these systems is in the hot rather than cold 
phase. This is consistent with the most X-ray luminous CGs being more mature, 
evolved systems. 

The relatively massive CGs (i.e., total stellar and Hi mass 
>10 n M q ) are typically more Hi-poor with 0.2% < /^ H i < 
6.6%, while the low-mass groups have 5.0% < < 40.6%. 

The most Hi-rich CG is HCG 30 with 40.7% of the baryonic 
mass in neutral hydrogen, while the most H i poor CG is HCG 
90 with <0.2% of the baryons in neutral gas, where the upper- 
limit is due to the non-negligible mass of hot gas that we have 
not included here. As the more massive CGs tend to be more 
X-ray luminous, it is possible the neutral gas in these systems 
was virialized to form the X-ray halo before the baryons could 
be used in star formation. 

4.5. Comparison of Diffuse X-Rays with Star Formation 

Desjardins et al. (2013) find that CGs with low specific star 
formation rates (sSFRs; i.e., the SFRs normalized to the galaxy 
stellar masses) are X-ray bright compared to CGs with higher 
sSFRs. This separates X-ray detected CGs into two types: (1) 
those with gas temperatures and luminosities consistent with 
virialization; and (2) those with hot gas associated with vigorous 
star formation (sSFR > 1.5 x 10 -11 yr -1 ). This latter class of 
star-forming groups appears to be a distinct class in L x - sSFR 
space. Diffuse emission is not detected in HCGs 7 and 22, which 
are groups with intermediate sSFRs. 

Using the UV+24 /xm SFRs from L. Lenkic et al. (2014, in 
preparation), which expands upon previous work by Tzanavaris 
et al. (2010), for the Walker et al. (2012) Expanded Sample of 
CGs, Figure 7 shows the diffuse X-ray luminosity as a function 
of the total group SFR. Note that SFR data are not available 
for HCGs 30, 40, and 68. Furthermore, the UV photometry of 
all of the galaxies in HCG 62 as well as galaxies HCGs 100C 
and D was unavailable, therefore we consider only IR SFR 
measurements for these galaxies. The purpose of combining the 
UV and 24 /xm SFRs is to ensure a complete census of star 
formation over the past ~100 Myr in both areas with low dust 
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Figure 7. CG diffuse X-ray luminosity as a function of the total group star 
formation rate as measured from the UV and 24 gm luminosities (Lenkic et al. 
2014, in preparation). There are no SFR data for HCGs 30, 40, and 68, and 
we do not include UV SFRs for HCGs 100C and 100D as these data were not 
available. An upper limit on the total SFR occurs when none of the galaxies in a 
group has detectable star formation. Groups with L x < 10 41 erg s -1 (HCGs 16, 
31, 59, and 90) tend to have galaxy-linked X-ray emission that increases with 
total group SFR, while the more X-ray luminous CGs have much lower total 
SFRs for their X-ray luminosity. 


column densities and star-forming regions enshrouded in dust. 
The lack of either a UV or IR SFR for a single galaxy is not 
expected to affect the SFR of that galaxy by a factor of >2 
for normal galaxies (excluding, e.g., LIRGs). We compare the 
diffuse X-ray emission against the SFRs rather than the sSFRs 
because it is the absolute rate of star formation that is most 
likely responsible for the X-ray emission when it is linked to 
individual galaxies. 

We find that the CGs in the L^-SFR space show two distinct 
distributions. At L x < 10 41 erg s -1 (HCGs 16, 31, 59, and 90), 
the hot gas is galaxy-linked, and the L x - SFR relation is found 
to be generally positive. This is expected if star formation is the 
source of the hot gas — more vigorous star formation leads to 
increased gas heating. Such heating may occur through various 
processes such as intense ionizing radiation from massive OB 
associations or supernova shocks. Further study of more groups 
with high total SFRs and low X-ray luminosities is required to 
better quantify this result. We note that the X-ray emission 
in HCG 90 is more likely due to tidal heating rather than 
star formation (Desjardins et al. 2013). For CGs with L x > 
10 41 erg s -1 , the SFR is generally lower in comparison. This 
overall decrease in SFR for X-ray bright systems may be 
attributed to gas stripping caused by the hot IGM (i.e., ram- 
pressure stripping; Gunn & Gott 1972), exhaustion of galaxy 
gas supplies to produce the hot IGM, or the cool gas supply was 
exhausted by star formation and is not necessarily coupled to 
the intragroup X-ray emitting gas. 

To investigate the star formation history, and how that relates 
to the X-ray luminosity of the groups, we also examine the 
g' — r', g' — i', and r' — i' optical colors of the first-ranked 
galaxy using data from Walker et al. (2013). We note that HCGs 
30, 42, 62, 68, and 90 are missing SDSS photometry. From the 
CGs with SDSS coverage, there are no obvious distinctions in 
color-color space between the first-ranked galaxies hosted in 
either X-ray luminous or X-ray non-detected CGs. There were 


also no correlations between the X-ray luminosity and any of 
the optical colors themselves. There is no evidence that the 
X-ray halos of the non-star-forming but X-ray luminous CGs 
were built-up by a recent (<1 Gyr) burst of star- formation in 
the first-ranked galaxy. 

5. SUMMARY AND CONCLUSIONS 

From a sample of 10 CGs, we detect five in X-rays using 
Chandra ACIS observations. We combine the CG X-ray tem- 
peratures and luminosities with those from Desjardins et al. 
(2013) to create a larger sample to examine the buildup of hot 
gas in CGs. The X-ray detected CGs in the combined sample 
range from 0.29-1.36 keV and to 4008-42 - 73 erg s -1 in temper- 
ature and X-ray luminosity, respectively. We then compared 
the X-ray properties against the stellar and Hi masses, line- 
of- sight velocity dispersions, and SFRs. Our results can be 
summarized as: 

1. From the cluster scaling laws (Section 4.2; Figure 3), we 
confirm our previous finding in Desjardins et al. (2013) 
that cool CGs with low velocity dispersions do not have 
cluster-like diffuse X-ray emission. However, we now find 
evidence that relatively massive and X-ray luminous CGs 
represent a population of groups that are consistent with 
the X-ray cluster scaling relations. This is consistent with 
previous X-ray studies of groups such as Connelly et al. 
(2012) and Lovisari & Reiprich (2013). 

2. In Section 4.3, we create a physically motivated definition 
of fossil groups based on the stellar mass distribution and 
the morphology of the most massive group galaxy, with 
no additional requirement for an X-ray bright halo. The 
stellar mass requirement follows from the description of 
a fossil group as an evolved system in which most of the 
mass is concentrated in one galaxy, which is supported by 
the work of Harrison et al. (2012). We require that: (1) 
the most massive galaxy contain >60% of the total group 
stellar mass; or, alternatively, the most massive galaxy 
contain >3 times more stellar mass than the next most 
massive galaxy; and (2) the most massive galaxy have an 
E/S0 morphology. This definition is effective at identifying 
low-mass fossil systems that may have been excluded by 
previous definitions (e.g., Jones et al. 2003). 

3. Using our fossil group definition, we identify 21 CGs 
that meet the stellar mass distribution criteria (Figure 5), 
but only six of these (HCGs 19, 22, 40, and 42, and 
RSCGs 44, and 86) host E/S0 galaxies as the most massive 
member, and are therefore fossil groups. The high fraction 
of groups (46%) with more than 60% of their stellar mass 
concentrated in the first-ranked galaxy is similar to the 
results of Connelly et al. (2012) for optically and X-ray 
selected groups. From the three CGs that have Chandra 
observations, only HCG 42 is X-ray luminous, and therefore 
may merge to form a massive, X-ray bright field elliptical. 
Conversely, HCGs 22 and 40 may merge to be X-ray faint 
field ellipticals. Furthermore, the spiral dominated groups 
that meet the stellar mass criteria for fossil groups (HCGs 
2, 4, 25, 26, 47, 56, 71, 79, and 100, and RSCGs 31, 34, 
64, and 66) 13 may represent interesting systems for study 
because they may be the dynamically young precursors to 
fossil systems. 


13 This list excludes HCG 31 because the aperture for the mid-IR photometry 
of 31 A includes 3 1C, and HCG 54 because this system is a false group and is 
instead made up of knots in a single galaxy. 
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4. From Sections 4.4 and 4.5, the X-ray luminosity of the CGs 
in our sample is correlated with the total stellar and H i mass 
of the groups (Figure 4), while the relationship with SFR 
is more complicated. At low values of L x , we qualitatively 
observe a positive trend between L x and SFR, while the 
opposite is true in X-ray bright systems (see Figure 7). Thus, 
X-ray luminous groups, especially those consistent with 
the cluster scaling relations, likely represent a population 
of mature, evolved systems. Low-mass groups are not a 
homogeneous population, and contain both dynamically 
young (e.g., HCG 31) and dynamically evolved (e.g., HCG 
22) systems. In general, low-mass groups, including those 
classified as fossil systems, are not X-ray luminous and fall 
well below the X-ray luminosities predicted by the cluster 
scaling relations. 

This last point raises the question of which mechanism takes 
precedence in the consumption of neutral gas in massive groups: 
star formation or IGM virialization. Specifically, when does 
the X-ray halo form in the evolutionary history of the group? 
The existence of HCGs 7 and 40, which represent relatively 
massive CGs with no evident X-ray emission and moderate 
SFRs, may indicate that the formation of the X-ray IGM 
occurs at a later evolutionary stage, supporting the hypothesis of 
Rasmussen et al. (2006) that X-ray underluminous loose groups 
are relatively dynamically young. HCG 16, another example of a 
high-mass group, and one that appears to be dynamically young 
(Konstantopoulos et al. 2013), was found to be dominated by 
galaxy-linked X-ray emission; however, a recent deep Chandra 
observation shows signs that there may be very faint diffuse 
X-ray emission in the intragroup space of HCG 16 (J. Vrtilek 
2013, private communication), and we refer the reader to the 
analysis of these data (E. O’Sullivan et al., in preparation). We 
note that in the case of HCG 7, the group has a low fraction of 
its mass in cool gas, therefore it may consume its gas reservoir 
prior in star formation before it is able to heat it to X-ray 
temperatures. Additionally, the development of a hot IGM could 
explain the anti-correlation of L x with SFR as the retardation 
in star formation due to ram-pressure stripping of cool gas in 
the group members. This puts our previous question into another 
perspective: when does star formation in galaxy groups end? The 
example of RSCG 3 1 , a group with a low SFR, low H i to stellar 
mass ratio, and high stellar to dynamical mass fraction, may be 
representative of a system in which star formation consumed 
the baryons and prevented the formation of an X-ray halo. 

Our findings indicate that potentially many low-mass galaxy 
groups are not X-ray luminous, and therefore X-ray surveys used 
to identify groups may create heavily biased samples that miss 
non-negligible fractions of the baryonic mass contained within 
groups of galaxies (see, e.g., Rasmussen et al. 2006; Connelly 
et al. 2012). Indeed, if the baryons are not in the molecular, 
Hi, or X-ray gas, then the evidence indicates that low-mass 
groups must contain a substantial fraction of baryons in either 
stars or a warm 10 5 -10 6 K gas phase. This warm gas would cool 
efficiently and be available again for star formation, but it would 
likely be replenished by gas heating in the group potential. We 
note that HCGs 30 and 31 have >30% of their mass in Hi; 
however, the next most Hi-rich CG contains only 11.3% of its 
baryonic mass in neutral gas, indicating that CGs with such 
an abundance of H i are exceptional rather than commonplace. 
Furthermore, it is unlikely that the stellar baryon fraction in 
low-mass groups is higher than in high-mass systems, as other 
studies find that the ratio of stellar to dynamical mass in groups 
is ~1% (Andreon 2010; Balogh et al. 2011; Connelly et al. 


2012). Thus, X-ray surveys only select the X-ray luminous, and 
therefore massive, groups while missing a significant population 
of low-mass groups. 
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